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Abstract. In this paper we analyze the quasi-stationary workload of a Levy-driven storage system. 
More precisely, assuming the system is in stationarity, we study its behavior conditional on the event that 
the busy period T in which time is contained has not ended before time t, as t — !> oo. We do so by 
first identifying the double Laplace transform associated with the workloads at time and time t, on the 
event {T > t}. This transform can be explicitly computed for the case of spectrally one-sided jumps. 
Then asymptotic techniques for Laplace inversion are relied upon to find the corresponding behavior in 
the limiting regime that t ^ oo. Several examples are treated; for instance in the case of Brownian input, 
we conclude that the workload distribution at time and t are both Erlang(2). 

Keywords. Levy processes * storage systems * quasi-stationary distribution * Laplace transforms * 
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1. Introduction 

Consider a storage system with Levy input, i.e., the process {Q{t))t that evolves as a Levy process 
X{t) that is reflected at 0. In mathematical terms, this means that the workload in the storage system 
at time t is given by 

Q{t) ■.= X{t)-iniX{s), 

s<t 

where Q{0) = x for some initial workload a: > 0. Assuming KX{1) < 0, there exists a stationary 
distribution tt of Q{t); it is seen that the stationary workload is distributed as the all-time supremum: 

7r(x) = P (supX(t) < x] . 
\t>o J 

In the sequel we add the subscript tt to the probability measure P and the associated expectation E 
when we wish to indicate that Q(0) is distributed according to this stationary distribution. 
Now let T denote the busy period, that is T = inf{t > : Q{t) = 0}. In this paper the object of our 
interest concerns the existence and characterization of the joint conditional distribution 

lim P^(g(0) e Ax, Q{t) edy\T>t)^: fi{dx, dy), 

t~¥OD 

where the convergence is to be understood in the weak sense. We study this so-called quasi-stationary 
distribution fi in detail; special attention is paid to the marginal distributions x R) = and 
/i(M X •) = /i^^(-). As an aside we mention that sometimes /i(K+ x dy) is called quasi-stationary 
distribution; here we do not follow that convention. 

A substantial body of work has been devoted to the analysis of quasi-stationary distributions. Over 
the past decades, various settings were considered; we here give a brief (non-exhaustive) overview. 
Seneta and Vere-Jones [TT], Tweedie IITSl , Jacka and Roberts |l9| focus on a Markov chain setting, Igle- 
hart |j8J addresses a random walk setup, Kyprianou IflOl considers the M/G/1 queue (i.e., a storage 
system with compound Poisson input), whereas Martinez and San Martin [ |13j treat the case of Brow- 
nian motion with drift. We also mention the contribution by Kyprianou and Palmowski ||T21 , who 
found the quasi-stationary distribution associated with a general light-tailed Levy process. Recently, 
Rivero fl6l (after appropriate scaling) found the quasi-stationary distribution for the specific situation 
that the Levy process under study has a jump measure with a regularly varying tail. 

The contribution of this paper is twofold. In the first place, a general formula for the double Laplace 
transform {Q{0),Q{t)) on the event {T > t}, that is, 

(1) / e-^X[e-"'?^°)-^'^(*\r>t]dt 

Jo 

is given (to which we refer to as the master formula). The derivation is based on the Wiener-Hopf 
factorization, and can be evaluated explicitly when all jumps are either all positive (the so-called 
spectrally-positive case) or all jumps are negative (spectrally-negative case). These formulae allow us 
to identify the quasi-stationary measures for the spectrally one-sided cases (relying on Tauberian-type 
theorems), which can be regarded as the second major contribution. 

The paper provides interesting insights into the distribution of the workload conditional on a long 
busy period. The distributions found tend to be stochastically larger than the normal, stationary dis- 
tribution. For instance in the case of regulated standard Brownian motion (with drift —1), both /i^^ (•) 
and Mi^^(') correspond to Erlang(2) distributions (with mean 2), whereas the stationary distribution 
is exponential (with mean i). This type of insights can potentially be used when setting up efficient 
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importance sampling algorithms ||6) (as in those algorithms a change of measure is looked for that 
mimics the distribution conditional on the rare event under consideration). 

This paper is organized as follows. In Section |2] preliminaries are given: (i) we first recapitulate a 
set of main results on fluctuation theory for Levy processes, and (ii) then present Tauberian theorems 
that are useful in the context of this paper (which can be used to identify the tail behavior of a random 
variable from its Laplace transform). The main objective of Section |3] concerns the derivation of the 
master formula, i.e., an expression for ([1) in terms of the Wiener-Hopf factorization (with explicit 
results for the spectrally one-sided cases). Our findings on the quasi-stationary distribution are then 
given in SectionH] The last section treats a number of examples. 



2. Preliminaries 

2.1. Levy processes. Here we follow flTl for definitions, notations and basic facts on Levy pro- 
cesses. Let in the sequel X = {X{t))t be a Levy process which is defined on the filtered space 
(fJ, {^t}t>o, P) with the natural filtration that satisfies the usual assumptions of right continu- 
ity and completion. Later if we write fx, it means that ¥x{X{Q) = a;) = 1 and Pq = P; similarly, E^; is 
expectation with respect to P^r. We denote by n(-) jump measure of X. Later we will focus on asym- 
metric Levy processes, which are either spectrally negative (having nonpositive jumps) or spectrally 
positive Levy processes (having nonnegative jumps). 

First passage times. For any Levy process we can define its Laplace exponent by 

(2) Ee''-^**) = e*'''(''\ 

for rj € Q such that the left hand side of (|2)l is well-defined (from now on we will assume that that this 
set 9 is not empty). Later, we also need the first passage time 

t{x) := niin{t : X{t) > x}, 

which for a spectrally negative process X with positive drift (i.e., KX{1) > 0) has Laplace transform 

(3) Ee~"^(^' = 6"*^")"^, 

where $(s) := sup{77 > : ?/'(7]) = s} is the right inverse of ■0 (see IITTl for details). 

Exponential change of measure. For e 9 we define a new probability measure P^ by the relation 

dPx 

we say that we have performed an exponential change of measure. Under P'', the process X is still a 
Levy process, but now with Laplace exponent 

(4) V'r,(/3) + 

We will use the subscript d to indicate that the quantity under consideration relates to P''. 

Dual process. We will also consider the so-called dual process Xt ~ ~Xt with jump measure n (0, y) — 
n {~y, 0). Characteristics of X will be indicated by using the same symbols as for X, but with a 
added. 
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Ladder heights. For the process X we define the associated {L~^(t), H(t))t: 

L-\t) 



inf{s > : L{s) > t} i{t< L(oo), 
oo otherwise. 



and 



H{t):=\ ift<i(oo), 
I oo otherwise, 



where L = {L{t))t is the local time at the maximum fill p. 140]. Recall that {L ^,H) is a bivariate 
subordinator with the Laplace exponent 



and with the jump measure Hh- In addition to this, we define the descending ladder height process 
{L^^{t), H{t))t>o with the Laplace exponent /3) constructed from the dual process X. Recall that 
under the stability assumption EX{1) < 0, the random variable L{oo) has an exponential distribution 
with parameter k(0, 0). Moreover, for a spectrally negative Levy process the Wiener-Hopf factoriza- 
tion states that 

(5) A^(^,/3)=<l>(^)+/3, k{^,(3)= l~f'^l 

- P 

see IHI p. 169-170]. It follows that k(0, 0) = i''{0+). 
We introduce a potential measure defined by 

/•oo 

'^(d.T,ds) = / F{L-\t) e ds,H{t) e dx) dt 
Jt=o 

with the Laplace transform /j^ e^''^"'^^^'^ {dx, ds) = 1/k((^, j3) and renewal function 

V{dx) = ^(da;, d.s) = E l{Hit)ed.}dt^ ■ 

In particular, 

(6) / e-^'^Vix) dx 



/3k(0,/3)- 

For a spectrally negative Levy process, the upward ladder height process is just a linear drift, and 
hence the renewal measure corresponds to the Lebesgue measure: 

(7) V{dx) ^ dx. 

From ^ we have Ull p. 195] that 

/•oo 



2.2. Tauberian-type results. Consider a function / : K M such that f{z) ~ for ^{z) < 0. Let 

f{z) := Jg°° e~^^f{x) dx be its Laplace transform. Consider its singularities; among these, let i^"^ < 
the one with the largest real part. Notice that this yields the integrability of |/(x)| da;. The 
inversion formula the reads 

-1 pa+ioo 
^ J a—ioo 

for some (and then any) a > i!)*. 



Quasi-stationary workload in a Levy-driven storage system 



5 



We now focus on a class of theorems that infer the tail behavior of a function from its Laplace trans- 
form, commonly referred to as Tauberian theorems. Importantly, the behavior of the Laplace transform 
around the singularity ■d'* plays a crucial role here. The following heuristic principle given in [IJ is 
often relied upon. Suppose that for -d*, some constants K and C, and a non-integer s > 0, 

/(i?) =K- Ci-d - + o((i9 - ry), as i d*. 

Then 

r(-s) 

where r(,s) is the gamma function. Below we specify conditions under which this relation can be 
rigorously proven. Later in our paper we apply it for the specific case that s = 1/2; recall that 

r(-i/2) = -2^. 

A formal justification of the above relation can be found in Doetsch HI Theorem 37.1]. Following 
Miyazawa and Rolski IT4|| , we consider the following specific form. For this we first recall the concept 
of the !21J-contour with an half-angle of opening 7r/2 < < tt, as depicted on |4j Fig. 30, p. 240]; also, 
'^c^* (ip) is the region between the contour 211 and the line 5R(z) = 0. More precisely, 

•^aiS) = {zeC; 5R(z) < 0, z 7^ a, I arg(z - a)| < S}, 

where arg z is the principal part of the argument of the complex number z. 

In the following theorem, conditions are identified such that the above principle holds; we refer to 
this as the Heaviside's operational principle, or simply Heaviside principle. 

Theorem 1 (Heaviside principle). Suppose that for f : C ^ C and ^* < the following three conditions 
hold: 

(Al) /(•) is analytic in a region ^c* W)for some n/2 < tp < -k; 

(A2) f{z) as \z\ ^oo for z€ ^^W; 

(A3) for some constants K and C, and a non-integer s > 0, 

(8) f{z) = K i{,>o} - c{z - cy + o{{z - cy), 

where 'Sq* (tp) 3 z ^ (^'^. 

Then 

r(-s) 

as a; oo, ivhere K := /(C*) ifs>0. 

We now discuss when assumption (Al) is satisfied. To check that the Laplace transform /(•) is analytic 
in the region '^^(^* (ip), we can use the concept of semiexponentiality of / (see |7||). 

Definition 2 (Semiexponentiality). / is said to be semiexponential if for some < < 7r/2, there exists 
finite and strictly negative 7(1?), defined as the infimum of all such a such that 

\f{e'^r)\ < e"'' 

for all sufficiently large r; here —4><i}<4> and sup 7(1?) < 0. 

Relying on this concept, the following sufficient condition for (Al) applies. 
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Proposition 3. [7", Thm. 10. 9f] Suppose that / is semiexponential with ^{1!}) fulfilling the following condi- 
tions: (i) 7 = 7(0) < 0, (ii) 7(1?) > 7(0) in a neighborhood of ■& — 0, and (iii) it is smooth. Then (Al) is 
satisfied. 

3. Master formula 

The objective of this section is to derive a general formula for the double Laplace-Stieltjes transforms 



/•OO 

Jo 



and 



Let be an exponentially distributed random variable with parameter -d > Q, independent of the 
process X. Denote 

:= supX(s), X{t) inf X(s). 

s<t 

Following the idea behind |2l Th. VI.20], we can prove the following result. 
Theorem 4. For a, /3, > 0, 

Proof. It is elementary that 

pOD poo 

Jo 

Recall that [llli Thm. 6.16(i)]) under P we have that X_{e.d) and X(e^) — X_{e^) are independent; in 
addition X(ei3) — X_{e{)) ~d X{ei}). It thus follows that 

= le-^*-E[e-^(^(^-)-^(^^»]E[e-'5^(^-),X(e^) > -x] 



^e-f^'=¥.[e-^^^''^]t[e^^^^^\X{ei)) < x\ 



Following flT] Th. 6.16(ii)], we obtain that 

and hence ^ 

E[e/5^(«^);X(e^) <x]= f e^'F(X{e.o) € dz) 

which completes the proof. □ 
We now evaluate L{i3; a, f3) for the spectrally one-sided cases. 
Proposition 5. IfX is a spectrally positive Levy process, then 



(10) ma,P) 
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Proof. Note that 

P(X(e^) >x)= P(t(x) < e^) = Ee-''^^"^' = g-i-W^. 
Integration by parts yields 

la ^ ^ - - ' 

The Pollaczek-Khintchine formula Illl Eqn. (4.14), p. 101] states that 

V''(0+)s 



^ e^^ P(X(e.) e dz) . (eC^-^C))^ - l) 



(11) nis):= e— ^(dx) 
It now follows that 

m a, P) = ^^1^ e^^P(X(e,) e dz) ^(dx) 

- ^..(^,/3)/3-$(^)io ^ 
?/''(0+) k(i?,0) $(i?) / a + l>(i?) a + /3 \ 

d Kid, 13) i3-^d) l^^(a + $(,?)) ^(a + /3) j ' 

The Wiener-Hopf factorization fill Section 6.5.2] and (|5) complete the proof. □ 
A similar result can be derived for the spectrally negative case. 
Proposition 6. IfX is a spectrally negative Levy process, then 

(12) L(.;.,/5)='^(^)-"-*W "^(O) 



<^{§)+l3 a + /? + $(0)?9- V'(a + $(0))' 

Proof. Recall the well-kown fact that 7r(d2:) ~ $(0)e~*(''^^ Ax by ||3ll. Applying Thm.|4]and interchang- 
ing the order of integration, 

^ 1 ^.(^,0) $(0) jg.^,(,+^(o))X(e„) 

79K(7?,/3)a + /3 + $(0) 

^ 1 ^.(^,0) $(0) ^^^(c,+^(0))X(e..,) 

?9K(t?,/3)a + /J + $(0) 
This gives by Eqns. l|5]l and l|9ll, in conjunction with the fact that d = K(d, 0)k(i9, 0), 

r(9 R) = i ^^(^,o) m ^(0,0) 

^ '"'^^ i9K(t?,/3)a + /3 + $(0)K(i?,a + $(0)) 

k(7?,0) 

^ K{d, /?) a + /3 + 4>(0) a + $(0)) 
$(i9)-a-$(0) $(0) 1 

$(z9) +/? a + /3-f $(0)i?-i/;(a + $(0))' 

This completes the proof. □ 

4. Quasi-stationary distribution 

In this section we use the Laplace transforms given in JTOl l and (121 to identify the quasi-stationary 
distribution /i(d.T, dy) for the spectrally one-sided cases. 
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4.1. Spectrally positive Levy process. We impose the following additional assumptions: 
[SPl] There exists < such that 

• < oo for 79+ < d, 

• '4>{'d) attains its strictly negative minimum at d* < 0, where < "d* < (and hence ip'{'d*) = 
0). 

Denote (* := ip{-d*) < 0. Note that the function $ can be considered in the complex domain. It is 
clearly analytic for 3ff(i?) > (*. But, as it will turn out, more is required to obtain the quasi-stationary 
distribution. 

[SP2] One can extend analytically a, (3) into (ip) for some tt/2 < ^ < tt. 

Example 7. Since $(■!?) is the Laplace exponent of a subordinator (viz. a first passage time process), 
we have the following spectral representation: 

/•oo 

(13) $(i?)=d+i9+/ (l-e-''")n+(dx), 

Jo 

and /o°°(a; A 1)11+ (dx) < oo ITTl Exercise 2.11]. From its definition we see that must be a singular 
point of ^{d). Moreover, if there exists a density of 11+ which is of semixponential type, then from 
Prop. |3] it follows then that $ is analytic in '^('■{(f>) and assumption [SP2] is satisfied. In particular, 
assumption [SP2] is for example satisfied for 

n+(dx-) = e^'^'x" dxl{^>o}, 

for a > -2. Clearly, then = C* cos i}. 

Theorem S.IfX is a spectrally positive Levy process satisfying conditions [SP1-SP2], then 

/i(dx,dy) = g+ (-i;'{0+)t^^h-^^] e^''^y-^^xV.,.{y)7Tidx)dy l{.>o,,>o}, 

where Q+ := {j^ e^*W-^<.{z)Az)-^. 
Corollary 9. We have 

V {H^W ) 

and 

iff {Ay) = Q+e''*^l/_^*(2/)d2/l{,>o}. 



Before we prove Thm.|8l we first present a few facts. Let fc* := y 2/V'"(t?*). 
Lemma 10. Under [SP1-SP2], 

$(1?) = r + fc*(i9 - C*)^/^ + o((i9 - C*)^^^) 

as I? ; C- 

Proof. From a Taylor series expansion and the condition that ij}' id*) = 0, we have 
After some rearranging, it is obtained that 
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We now substitute t9 = ^{s), and use ilj{<t{s)) ^ s to complete the proof. □ 
Proposition 11. Under [SP1-SP2] we have 

lim E,[e-"«(°)-^'3W|T>t] 

Proof. By [SP1-SP2] and Prop.|5]L(i?, a, /3), as given in ((TOl l as a function of i?, is analytic in (0) for 
7r/2 < < TT when 'ip{a + ^(z)) is analytical there. Recall that is analytic in this region and note 
that %lj{a + $(2:)) is analytical there since = z is analytical in this region. Thus condition (Al) 

of Thm.[T]is satisfied. 

To check that condition (A2) of Thm.[T]holds for L{d, a, /3), it suffices to prove that 



(15) 



1 a + ^{z) 



tends to for z G Qc^* {(f>) tending to 00 two-dimensionally, that is in particular for 3 z ±00. From 
(131 it follows that b := ^<t{z) ±00. Similarly one can prove that either a := a + 5R$(z) < is 
constant or tends to —00. Then 



-0(a + 6i) _ 1 
\a + bi\ " + 

~ Vi'(O) - / xll{dx) + acr 



^-0'(O)& + 06(7^ + J (sinba: - 51{|j,|<i}a;)n(da:) 





2 



which is either bounded or tends to —00 when a is constant or a ^ —00. Taking into account the term 
1/z in i fTSl l, this completes the verification of condition (A2). We will check now that also condition 
(A3) of Thm.Uis satisfied. Now using Lemma [TOl we write 

= ij{a + + i''{a + r )k*(i9 - C)^^^ + o((t9 - CY^^). 
Hence by Prop. |5] (for | C*) we have, for some K, 

?/''(0+) / a + (3 a + i!)* + k*{^-C)^/^ + o{{d -CY/^)\ 



= K- 



K 



i9 - '/'(a + + i'' (a + ^*)k*{'d - O^/a + o{{§ - C^Y'^) 

V;'(0+) a + §* + k*{d - C*)^/^ + o((i? - C)i/2 

1} - V'(/3) V'(a + 1?^) + i/i'(a -f 79*)A:*(i9 - C)i/2 + o{{^ - C)^/^) 

ip{a + 'd*) - V)'(a + i9*)fc*(i?-C*)i/2 + o((t?-C*)i/2) ' 

Thus we obtain that, for some K, 

Conclude by invoking 'Heaviside' that 

E^f,-.g(o,-.g,)^^ ^ ^ ^-(0+)fc-^(a + ^-)-(a + ^-)^-(a + ^-) 0(i + ,(1)). 
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By setting a = (3 = we have 

iCy r(-i/2)- 

It is now seen that Eqn. ((Til l holds and the proof is completed. □ 



Proof of Theorem\8\ From Prop. [TT] it follows that 



OO /'OO 



where 



fi{a,P)= / e-"^e-'5V(da:,d2/)=i+(a)^+(/3), 
Jo Jo 

^+(")— C TT, ' ^+(P) 



• By the Pollaczek-Khintchine formula ((TTt and fill (4.14), p. 101] applied for the dual, we derive 
Hence, has the desired form. 

• The dual process X is spectrally negative, so that V-^* {y) = y and (0, 13) = ^ for all ?^ > 0. The 
Wiener-Hopf factorization gives (up to a multiplicative constant k that relates to the normalization of 
the local time) that under P^'' for all G R we have 

(16) V-i^* W = (0, -I?) 

for all i9 < — From (0) and ||6ll we have that 

e-^V*^F_^*(y) dy 



" " (/3-^*)a^_^*(0,/3-^9*) 

k k k 



Conclude that (') has the desired form, which completes the proof. □. 

Remark 12. The transform A+( ) can be used to interpret the quasi-stationary distributions. Because 
of PoUaczek-Khinchine, 

is a Laplace transform (i.e., corresponding to an exponentially twisted version of the steady-state 
workload). In addition, by virtue of |5l Lemma 3.5], 

2 a'ip'{a) — 'ip{a) 



is a Laplace transform, and therefore also its i?*-twisted version 



(17) 



This reasoning indicates that, conditional on a long busy period, Q(0) is distributed as the sum of 
three independent random variables. Two of these are distributed as the i9*-twisted version of the 
steady-state workload, while a third has transform ((TTl l. <) 
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4.2. Spectrally negative Levy process. Like for the spectrally positive case, also in the spectrally 
negative case we need to impose additional assumptions to find the quasi-stationary distribution. 
[SNl] There exists i?_ > such that 

• V(t?) < oo for < I? < "d^, 

• '4^{'d) attains its strictly negative minimum at d* > 0, where < -d* < (and hence tp'{'d*) = 
0). 

[SN2] $ is analytical in {(j)) for 7r/2 < < tt, where C* := ^{^*) < 0. 

Example 13. Since for spectrally negative Levy process is the Laplace exponent of a subordinator 
(viz. a first passage time process), the spectral representation 

/•oo 

(18) $(i9)=d_i9+/ (1 - e-''"^)n_(da;), 

Jo 

applies, with Jf^{x A l)n_ (dx) < oo; cf. ([131 . This means that if there exists a density of n_ which is 
of semixponential type, then Prop. |3] entails that $ is analytic in (0) and hence assumption [SN2] 
is satisfied. 

Theorem 14. If X is a spectrally negative Levy process satisfying conditions [SN1-SN2], then 

Ai(dx,dy) = Q_(r)2ye-''*(-+^)e-*W^K,.(x)da;d2/ l{,>o,,>o}, 
where := [f^ e~^'^^°'^+^*>-Vi,*{z) Az)-^. 
Corollary 15. We have 

/if (dy) = g_e-(*("'+''*'-K,.(x)dxl{,>o} 

and 

Observe that (•) corresponds with an Erlang(2) distribution. The proof of these results is based on 
the following lemma, which is proven as Lemma [TOl 

Lemma 16. Under [SN1-SN2], 



asd i C*/ ^here k* yy2/tp"{d*). 

Proof of Theorem [l4l Note that all assumptions of Thm.[T]are satisfied by Prop.lH In particular, as\z\ 
tends to infinity in {(j)) function | $(2;) | is either bounded or tends to infinity. In both cases condition 
(A2) is satisfied. Moreover, for some K, 

as i9 I C,* , and we can conclude by 'Heaviside' that 

(t^* + /?)2^(a + $(0))-Cr(-l/2)^ l^ + oi^)) 

as t — > 00. Therefore, 

/i(a,/3) = lim E^ie-^'SW-WW |y > ^ A_{a)B_(l3), 

where 

-C* ~ fi?*)^ 



V'(a + $(0))-C' (i9*+/3)2' 
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It is not hard to see that the proposed density indeed corresponds with this transform. □ 

5. Examples 

In this section we illustrate our theory by means of a number of examples. We indicate for what Levy 
processes our assumptions are fulfilled, and for a few of those processes we perform the computa- 
tions. 

According to Vigon's theory of philanthropy ||T91 , a (killed) subordinator is called a philanthropist if 
its Levy measure has a decreasing density on M.^ . Moreover, given any two subordinators Hi and H2 
which are philanthropists, providing that at least one of them is not killed, there exists a Levy process 
X such that Hi and H2 have the same law as the ascending and descending ladder height processes 
of X, respectively. Suppose we denote the killing rate, drift coefficient and Levy measures of Hi and 
H2 by the respective triples (6, 6, Uh^ ) and (6, S, Hh-^)- Then (I9l shows that the Levy measure of X 
satisfies the following identity 

/•oo 

n(x, 00) = / IIh2{u,oo)I1h-^{x + du) + 5nH-^{x) + bIlHi{x,oo), x > 0, 
Jo 

where tthi {x) is the density corresponding to . By symmetry, an obvious analogue of the above 
equation holds for the negative tail n(— 00, x), with a; < 0. 

Choosing then e.g. H2{t) = t and Hi{t) = T{t) with Laplace exponent $ and jump measure 11+ being 
semiexponential, then, using the above construction, we can easily give examples of a spectrally 
positive Levy process X satisfying conditions [SP1-SP2]. Similarly, using the above method we can 
construct spectrally negative Levy processes satisfying [SN1-SN2]. 

Usually these conditions can be verified in a straightforward manner, as we did in the examples 
below. 

Example 17. M/M/1 queue. In this case 

Nit) 

(19) X{t) = ^ <7, - t, 

1=1 

where ct^ (where i — 1,2, ...) are i.i.d. service times that have an exponential distribution with mean 
l/u. The arrival process is a homogeneous Poisson process N{t) with rate A; it is assumed that g :— 
Xjv < 1. We apply the theory of Section l4T] 
We have 



■4'iv) = - A 1 ' — = ?/ 

V V + J rj + I' 

yielding -d* = \f\v — v, and = —{\fv — VA)^- Furthermore 



= ^ 

and hence assumptions [SP1-SP2] are satisfied. 

Let us first concentrate on/ig^(-). We use Remark [I2I Using that 

we obtain 

(a + i?*)V;'(a + 77*)-V>(Q! + 79*) (7?^)2 / 
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This corresponds with the sum of two Exp(\/Ai^) random variables. Also, 



a + vAi^ 
a + \f\v — A 



a + V Aj^ 



corresponding with a shifted-Geom(y^)-distributed number of Exp(\/AZ') random variables. Con- 
clude that fJ.^^{-) corresponds to the sum of M independent Exp{\/Xi') random variables, where 

i.e., il/ has a negative binomial distribution with parameters 2 and A similar form is found for 
the general light-tailed M/G/1 case. 

Let us now study It is a matter of straightforward calculus to find that 

B+{P) = {v + P) 



A partial fraction expansion argument gives that this equals 



V - \J\v + /? 

In other words, the quasi-stationary distribution at time t (for t large) equals a mixture of an expo- 
nential and an Erlang(2) distribution. <0> 



Example 18. Linear Broionian motion. In this case X{t) = aB{t) — t, where cr > and B{t) is a standard 
Brownian motion. Remark that this process is spectrally positive and spectrally negative, so we can 
use both Thm. |8]and Thm. [ll 

Let us first see what the spectrally positive results would give. It is not hard to check that 

so that, in the setting of Section |4IT1 -d* = — l/cr^ and (* = —1/{2<t'^). It is a matter of straightforward 
computations now to obtain that 



Conclude that the quasi-stationary distributions of Q(0) and Q{t) (t large) are both Erlang(2) with 
mean whereas the stationary workload has an exponential distribution with mean 1/{2<t'^). (In 
the decomposition of Remark[12l the first two random variables have exponential distributions with 
mean I/ct^, the third is equal to 0). Interestingly, the relation with the Erlang(2) distribution has also 
been observed in, e.g., ISl lTSllTSl . 

The same result can be obtained by using the results from Section 14.21 Now i?* = l/cr^ and C* = 
-1/(2ct2). It is easily checked that $(0) = 2/cr2. As expected, we obtain Ji_ (a) = andi?_(/3) = 

B+W). 

In fact, in this case the quasi-stationarity distributions can be found in an explicit manner. For ease 
we restrict ourselves to studying just we do so by investigating the density 
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We rely on the standard equality 

and the fact that Q(0) (unconditioned) has an exponential distribution with mean i; here, $n(') de- 
notes the distribution function of a standard Normal random variable, where 5'N(a;) := 1 — ^N(a;) is 
its tail. It is known that, as x oo, 

(20) vi,^(,)^fi-J_+3\ 



27r 

Let us first determine PttCT > t) (where Q(0) has an exponential distribution with mean i), which 
can evidently be rewritten as 

Consider the first integral of the previous display. It can be evaluated as 

oo \/2n J -oo Jo v 27r 

Likewise, the second integral can be rewritten as 

1 1 2 /""^ f-t-vVt I ^ ^ 



J-oo v27r J-oo Jo v27r 

J-oo v27r V27r 

We arrive at 

(21) P,(r > = *N(\/i) + *n(-%/<) + 2i$N(-\/t) - 2\/< • -^e-^*. 

V 27r 

Using (|20] | it is readily verified that, for t large, 

¥^{T > t) ^ A . -^e-i*. 

Also, it holds that 

(22) ^¥^^Q^<q^T>t)^2er^''^^(^—/)-2<P 



so that we now have an explicit expression for ft {q), viz. the ratio of l|22t and (|2Tl l. Due to the asymp- 
totic equivalence (|20] |, Expression (|22|l behaves for t large as 



V2^ \t-q t + qj \ 2 t J \tVt V^^ 
We conclude that we again find that the quasi-stationary distribution of Q(0) is Erlang(2) with ex- 
pected value 2. 
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QUASI-STATIONARY WORKLOAD 
IN A LEVY-DRIVEN STORAGE SYSTEM 



MICHEL MANDJES, ZBIGNIEW PALMOWSKI, AND TOMASZ ROLSKI 



Abstract. In this paper we analyze the quasi-stationary workload of a Levy-driven storage system. 
More precisely, assuming the system is in stationarity, we study its behavior conditional on the event that 
the busy period T in which time is contained has not ended before time t, as t — !> oo. We do so by 
first identifying the double Laplace transform associated with the workloads at time and time t, on the 
event {T > t}. This transform can be explicitly computed for the case of spectrally one-sided jumps. 
Then asymptotic techniques for Laplace inversion are relied upon to find the corresponding behavior in 
the limiting regime that t ^ oo. Several examples are treated; for instance in the case of Brownian input, 
we conclude that the workload distribution at time and t are both Erlang(2). 

Keywords. Levy processes * storage systems * quasi-stationary distribution * Laplace transforms * 
Heaviside principle * fluctuation theory 
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1. Introduction 

Consider a storage system with Levy netput, i.e., the process {Q{t))t that evolves as a Levy process 
X{t) that is reflected at 0. In mathematical terms, this means that the workload in the storage system 
at time t is given by 

Q{t) ■.= x + X{t) - inf (x + X{s)y, 

s<t 

where — min{a, 0}. Assuming KX{1) < 0, there exists a stationary distribution tt of Q{t); it is 
known that the stationary workload is distributed as the all-time supremum: 

tt{x) =P (supX(t) < x] . 
\t>o J 

In the sequel we add the subscript tt to the probability measure P and the associated expectation E 
when we wish to indicate that (5(0) is distributed according to this stationary distribution. 
Now let T denote the busy period, that is T = inf{i > : Q{t) = Q]. In this paper the object of our 
interest concerns the existence and characterization of the joint conditional distribution 

lim P^(g(0) G dx, Q{t) edy\T>t)^: fi{dx, dy), 

where the convergence is to be understood in the weak sense. We study this so-called quasi-stationary 
distribution ji in detail; special attention is paid to the marginal distributions x M) = and 
^(M X •) ~ /i^^(-). As an aside we mention that sometimes x dy) is called quasi-stationary 

distribution; here we do not follow that convention. 

A substantial body of work has been devoted to the analysis of quasi-stationary distributions. Over 
the past decades, various settings were considered; we here give a brief (non-exhaustive) overview. 
Seneta and Vere-Jones ITTII , Tweedie IITSl , Jacka and Roberts ||9l focus on a Markov chain setting, Igle- 
hart |8| addresses a random walk setup, Kyprianou ITOl considers the M/G/1 queue (i.e., a storage 
system with compound Poisson input), whereas Martinez and San Martin IITSl treat the case of Brow- 
nian motion with drift. We also mention the contribution by Kyprianou and Palmowski 11121 , who 
foimd the quasi-stationary distribution associated with a general light-tailed Levy process. Recently, 
Rivero flSl (after appropriate scaling) foimd the quasi-stationary distribution for the specific situation 
that the Levy process under study has a jump measure with a regularly varying tail. 

The contribution of this paper is twofold. In the first place, a general formula for the double Laplace 
transform ((3(0), Q(t)) on the event {T > t}, that is, 

(1) / e-''X[e""^(°'-^'3(t)^T>i]df 

is given (to which we refer to as the master formula). The derivation is based on the Wiener-Hopf 
factorization, and can be evaluated explicitly when all jumps are either all positive (the so-called 
spectrally-positive case) or all jumps are negative (spectrally-negative case). These formulae allow us 
to identify the quasi-stationary measures for the spectrally one-sided cases (relying on Tauberian-type 
theorems), which can be regarded as the second major contribution. 

The paper provides interesting insights into the distribution of the workload conditional on a long 
busy period. The distributions found tend to be stochastically larger than the normal, stationary dis- 
tribution. For instance in the case of regulated standard Brownian motion (with drift —1), both y.^^{-) 
and Mi^^(0 correspond to Erlang(2) distributions (with mean 2), whereas the stationary distribution 
is exponential (with mean i). This type of insights can potentially be used when setting up efficient 
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importance sampling algorithms ||6) (as in those algorithms a change of measure is looked for that 
mimics the distribution conditional on the rare event under consideration). 

This paper is organized as follows. In Section |2] preliminaries are given: (i) we first recapitulate a 
set of main results on fluctuation theory for Levy processes, and (ii) then present Tauberian theorems 
that are useful in the context of this paper (which can be used to identify the tail behavior of a random 
variable from its Laplace transform). The main objective of Section |3] concerns the derivation of the 
master formula, i.e., an expression for ([1) in terms of the Wiener-Hopf factorization (with explicit 
results for the spectrally one-sided cases). Our findings on the quasi-stationary distribution are then 
given in SectionH] The last section treats a number of examples. 

2. Preliminaries 

2.1. Levy processes. Here we follow flTl for definitions, notations and basic facts on Levy pro- 
cesses. Let in the sequel X = {X{t))t be a Levy process which is defined on the filtered space 
(O, {^t}t>o, P) with the natural filtration that satisfies the usual assumptions of right continu- 
ity and completion. Later if we write P^, it means that Pj;(X(0) = x) = 1 and Pq = P; similarly, Ej; 
is expectation with respect to P^;. We denote by n( ) the jump measure of X. Later we will focus 
on asymmetric Levy processes, which are either spectrally negative (having nonpositive jumps) or 
spectrally positive Levy processes (having nonnegative jumps). 

First passage times. For any Levy process we can define its Laplace exponent by 

(2) Ee''^^*) = e*'''^''), 

for ?7 G 9 such that the left hand side of l|2]l is well-defined. Later, we also need the first passage time 

t{x) := min{t : X{t) > x}, 
which for a spectrally negative process X with positive drift (i.e., KX{1) > 0) has Laplace transform 

(3) Ee-"^(^) = e-*(")=^, 

where $(5) mi{ri > : ip{ri) > s} is the right inverse of ijj (see flT] for details). 

Exponential change of measure. For e 9 we define a new probability measure P^ by the relation 

^ gr,(X(t)-:c)-^(T,)t . 

dP. 

we say that we have performed an exponential change of measure. Under P'', the process X is still a 
Levy process, but now with Laplace exponent 

(4) V^(/3) + 

We will use the subscript 77 to indicate that the quantity under consideration relates to P**. 

Dual process. We will also consider the so-called dual process Xt = —Xt with jump measure n (0, y) = 
n {—y, 0). Characteristics of X will be indicated by using the same symbols as for X, but with a 
added. 

Ladder heights. For the process X we define the associated {L^^{t),H{t))t: 

j^-i^^-^ ._ I inf{s > : L{s) > t] if i < L(oo), 
I 00 otherwise. 
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and 

H{t) := I ift<i(oo), 
I oo otherwise, 

where i = {L{t))t is the local time at the maximum flTl p. 140]. Recall that {L^^,H) is a bivariate 
subordinator with the Laplace exponent 



and with the jump measure Il^f . In addition to this, we define the descending ladder height process 
{L^^{t), H{t))t:^Q with the Laplace exponent k{'d, f3) constructed from the dual process X. Recall that 
under the stability assumption EX{1) < 0, the random variable L(oo) has an exponential distribution 
with parameter k(0, 0). Moreover, for a spectrally negative Levy process the Wiener-Hopf factoriza- 
tion implies that 

(5) = 

see HH p. 169-170]. It follows that k(0, 0) = ^'{0+). 
We introduce a potential measure ^ defined by 

/•oo 

'^(da;,ds)= / P(L"i(t) G ds,i7(t) e dx) dt 
with the Laplace transform Jjp e~'^'^~^^'^{dx, ds) = 1/k{(p, /?) and renewal function 

V{dx) = ^ ^ ^(dx, ds) = IE (^^ ^ l{Hmd.}dt^ ■ 

In particular, 

(6) / e-^^F(x) dx 



/3k(0,/3)- 

For a spectrally negative Levy process, the upward ladder height process is just a linear drift, and 
hence the renewal measure corresponds to the Lebesgue measure: 

(7) V{dx) ^ dx. 

From ^ we have Ull p. 195] that 

/•oo 

According to our convention K?, F^, k^, are quantities computed under P''. 

2.2. Tauberian-type results. Consider a function / : M R such that /(x) = for x < 0. Let f{z) := 

Jp°° e~^^/(x) d.x for z G R be its Laplace transform. Consider singularities of f{z); among these, let 
i3* < the one with the largest real part. Notice that this yields the integrability of \f{^) \ dx. The 
inversion formula reads 

-I pa+ioo 

/(^) = — / /»e^"dz 

^ J a— ioo 

for some (and then any) a> 

We now focus on a class of theorems that infer the tail behavior of a function from its Laplace trans- 
form, commonly referred to as Tauberian theorems. Importantly, the behavior of the Laplace transform 
around the singularity i?* plays a crucial role here. The following heuristic principle given in f\\ is 
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often relied upon. Suppose that for d*, some constants K and C, and a non-integer s > 0, 

/(i?) = A' - C(i? - i)*Y + o((i9 - i}*y), as I? ; 

Then 

r(-s) 

where r(,s) is the gamma function. Below we specify conditions under which this relation can be 
rigorously proven. Later in our paper we apply it for the specific case that s = 1/2; recall that 

r(-i/2) = -2^. 

A formal justification of the above relation can be found in Doetsch HI Theorem 37.1]. Following 
Miyazawa and Rolski fT4ll , we consider the following specific form. For this we first recall the concept 
of the SlU-contour with an half -angle of opening 7r/2 < ^ < tt, as depicted on [Aj, Fig. 30, p. 240]; also, 
is the region between the contour 211 and the line 3? (z) = 0. More precisely, 

^„(V') = {z e C;3fi(z) < 0,Z7^a,|arg(z-a)| < V}, 

where arg z is the principal part of the argument of the complex number z. In the following theo- 
rem, conditions are identified such that the above principle holds; we refer to this as the Heaviside's 
operational principle, or simply Heaviside principle. 

Theorem 1 (Heaviside principle). Suppose that for / : C — )- C and ^* < the following three conditions 
hold: 

(Al) /(•) is analytic in a region {"i') for some 7r/2 < tp <t:; 

(A2) f{z) -^Oas \z\ oo for ze ^^W; 

(A3) for some constants K and C, and a non-integer s > 0, 

(8) f{z) = K l{,>o} - C{z - Cr + o{{z - CY). 

where 'Sq* (tp) 9 z — > C*- 

Then 

as X ^ CO, where K := /(C*) ifs>0. 

We now discuss when assumption (Al) is satisfied. To check that the Laplace transform /(•) is analytic 
in the region '^^(^* (ip), we can use the concept of semiexponentiality of / (see iTl). 

Definition 2 (Semiexponentiality). / is said to be semiexponential if for some < < n/2, there exists 
finite and strictly negative 7(7?), defined as the infimum of all such a such that 

\f{e''^r)\ < e'"' 

for all sufficiently large r; here —(p <'d<(p and sup7(?9) < 0. 

Relying on this concept, the following sufficient condition for (Al) applies. 

Proposition 3. [7, Thm. 10.9f] Suppose that / is semiexponential with ^{d) fulfilling the following condi- 
tions: (i) 7 = 7(0) < 0, (ii) 7(??) > 7(0) in a neighborhood of-d = 0, and (iii) it is smooth. Then (Al) is 
satisfied. 
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3. Master formula 

The objective of this section is to derive a general formula for the double Laplace-Stieltjes transforms 

/>oo 

L,(z^;a,/3) := / e-''%[e-"^-^0(*', T > t] 
Jo 

and 

/■oo />oo 

L{§-a,l3):= a, /3) dP(Q(0) < x) = / e-''*E^[e-"^(°)-'30(*), T > i] di. 

Let e§ be an exponentially distributed random variable with parameter -d > Q, independent of the 
process X. Denote 

:= supX(s), X{t) := inf X(s). 
From IIT2I eq. (10)] we can get the the following result. 
Theorem 4. For a, /3, > and x>Q, 

We now evaluate L{d] a, (3) for the spectrally one-sided cases. 
Proposition 5. IfX is a spectrally positive Levy process, then 

Vi'(0+) / a + P a + $(i9) A 



(9) L{^;a,l3) 

d - ij{/3) \ij{a + /3) + 

Proof. Note that 

P(X(e^) >x)= V{t{x) < e^) = Ee'^^^^^ = e^^^*)" 
Integration by parts yields 

'e''^P(X(e.)edz) = -i§-(e(^-*(''))^-l) 

The Pollaczek-Khintchine formula Illl Eqn. (4.14), p. 101] states that 

Vi'(0+)s 



(10) Hs) / ^(dx) 

Jo ij(s) 

It now follows that 

m /?) = e"^"^^'" £ e^^P(X(e,) € dz) nidx) 

i?K(z?,/3)/3_$(t?) io ^ y ^ ^ 

Vi'(0+) k(i?,0) $(i9) / a + $(i?) a + ^ \ 

1? k(i?,/3) l^V>(a + $(i9)) V'(a + /3) / ' 

The Wiener-Hopf factorization (llj Section 6.5.2] and (|5) complete the proof. □ 

A similar result can be derived for the spectrally negative case. 

Proposition 6. IfX is a spectrally negative Levy process, then 

$(7?)-a-$(0) $(0) 1 



(11) ma,p)^ 



$(i9)+/3 a + ;3 + $(0)i9-?/'(a + $(0))' 
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Proof. Recall the well-kown fact that 7r(da;) = $(0)e by ||3ll- Applying Thm.|4]and interchang- 

ing the order of integration, 

^ 1 k(^?,0) $(o) 

i?K(i?,/3)a + /3 + $(0) 

^ 1 .^(7?, 0) $(0) ^^_(^+^(o))x(e,) 

This gives by Eqns. ^ and fill Th. 6.16(ii)], in conjunction with the fact that i? ~ K{d, O)k{'0, 0), 

1k(79,0) $(0) k(0,0) 



(3) 



1? K(t?, /3) a -f /3 -f $(0) a + $(0)) 
$(i9)-a-$(0) $(0) 1 



^{'d)+l3 a + /3 + $(0)i9-'(/'(a + *(0))' 
This completes the proof. □ 



4. Quasi-stationary distribution 

In this section we use the Laplace transforms given in ||9) and l(TT| l to identify the quasi-stationary 
distribution fi{dx, dy) for the spectrally one-sided cases. 

4.1. Spectrally positive Levy process. We impose the following additional assumptions: 
[SPl] There exists < such that 

• ijj{-d) < oo for i9+ < I?, 

• -0(1?) attains its strictly negative minimum at < 0, where ?9+ < t?* < (and hence = 
0). 

Denote C* := < 0. Note that the function $ can be considered in the complex domain. It is 

clearly analytic for 3fi(i?) > (*. But, as it will turn out, more is required to obtain the quasi-stationary 
distribution. 

[SP2] One can extend analytically L{-d; a, (3) into (ip) for some 7r/2 < ^ < tt. 

Example 7. Since <l>(i?) is the Laplace exponent of a subordinator (viz. a first passage time process), 
we have the following spectral representation: 

/•oo 

(12) m)^d+^+ (l-e-''")n+(dx), 

Jo 

and /q°°(x a l)n+(da;) < oo [Til Exercise 2.11]. From its definition we see that (* must be a singular 
point of $(i?). Moreover, if there exists a density of 11+ which is of semiexponential type, then from 
Prop. |3] it follows then that $ is analytic in ^^(j*((/)) and assumption [SP2] is satisfied. In particular, 
assumption [SP2] is for example satisfied for 

(13) n+(dx-) = e«*^a;"da;l{,>o}, 

for a > —2. Clearly, then 7(1?) = (* cosi?. Note that ([131 is satisfied for a linear Brownian motion 
X{t) = B{t) — t. Indeed, in this case ^{-d) is the Laplace exponent of the inverse Gaussian process 
with Levy measure n+(da;) = {2'7t)~^^'^x~^^'^ exp{— |} dx. 
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Theorem 8. IfX is a spectrally positive Levy process satisfying conditions [SP1-SP2], then 



where Q+ {j^ e^*''V-^*{z)dz)-^. 
Corollary 9. We have 

and 

fifidy) = Q+e''*yV.^.{y)dyl{y>o}- 



Before we prove Thm.|8l we first present a few facts. Let k* := y 2/-0"(i?*). 
Lemma 10. Under [SP1-SP2], 

= I?* + k*{^ - c*)i/^ + o((i? - cy^^) 

asdi C- 

Proof. From a Taylor series expansion and the condition that ?/;'(t?*) — 0, we have 
After some rearranging, it is obtained that 



We now substitute d ~ and use ?/'(<l>(s)) = s to complete the proof. □ 
Proposition 11. Under [SP1-SP2] we have 

lim E,[e-"«(°)-^«W|T><] 

V ^P^ia + ^*) J \C- 4,(13)) ■ 

Proof. By [SP1-SP2] and Prop.|5]L(i?, a, (3), as given in (|9]l as a function of i?, is analytic in ^^^* (0) for 
7r/2 < (f) < TT when ?/'(q! + $(z)) is analytical there. Recall that $ is analytic in this region and note 
that ■(/;(« + $(z)) is analytical there since t/j{'t{z)) = z is analytical in this region. Thus condition (Al) 
of Thm.[T]is satisfied. 

To check that condition (A2) of Thm.[T]holds for L{-d, a, (3), it suffices to prove that 



(15) 
tend 

ing ^{z) := a + bi note that 



1 a + ^{z) 



z tp{a + $(z)) 

tends to for z e CJ^. (</>) tending to 00 two-dimensionally, that is in particular for 3 z — ±00. Denot- 



Vi(a + h\) _ 1 



('4,'{0)b + ab(j'^+ I (sin6x - 61{|^|<i}a::)n(da::) 
\ J —00 



\a + bi\ + 52 

cannot be zero. Taking into account the term 1 /z in |(T5] |, this completes the verification of condition 
(A2). We will check now that also condition (A3) of Thm. [T]is satisfied. Now using Lemma [TOl we 
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write 

= i>{a + r) + f (a + r)fc*(i? - Cf'' + o((i? - cf''). 
Hence by Prop. |5] (for d !^ C,*) we have, for some K, 

?A'(0+) / a + P a + d* + k*(§^C*f/'^ + o((^-CY/'^) 



K 



1? - ?/'(^) tP{a + I?*) + + d*)k*{d - 01/2 + o((?9 - C)^/^) 

Tpia + -i^'{a + 'd*)k*{^ - c*y/^ + o{{i9 - c*y/^) ' 

Thus we obtain that, for some K, 
Conclude by invoking 'Heaviside' that 

^ ^ iP{a + i)*) r(-i/2) ^ ^ 

By setting a ~ p = {] we have 

It is now seen that Eqn. |[T4t holds and the proof is completed. □ 
Proof of Theorem^ From Prop. [TT] it follows that 

noo 

where 



A+{a):=C r— — , B+(l3) — 



By the Pollaczek-Khintchine formula dlOl l and HTl (4.14), p. 101] applied for the dual, we derive 



00 



Hence, a^2^(-) has the desired form. 

• The dual process X is spectrally negative, so that V^^* (y) ~ y and (0, ■&) — ■& for all ?9 > 0. The 
Wiener-Hopf factorization gives (up to a multiplicative constant k that relates to the normalization of 
the local time) that under P^'' for all G R we have 

(16) ij-i,* (i?) -k§K_i). (0, -I?) 
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for all ?? < — From lHJl and l|6| we have that 



oo -1 



k k k 



Conclude that (• ) has the desired form, which completes the proof. □. 

Remark 12. The transform j4+( ) can be used to interpret the quasi-stationary distributions. Because 
of Pollaczek-Khinchine, 

is a Laplace transform (i.e., corresponding to an exponentially twisted version of the steady-state 
workload). In addition, by virtue of O Lemma 3.5], 

2 a?/)' (a) - ■0(a) 



is a Laplace transform, and therefore also its i?*-twisted version 



(17) 



This reasoning indicates that, conditional on a long busy period, Q(0) is distributed as the sum of 
three independent random variables. Two of these are distributed as the i9*-twisted version of the 
steady-state workload, while a third has transform (flTt . (} 

4.2. Spectrally negative Levy process. Like for the spectrally positive case, also in the spectrally 
negative case we need to impose additional assumptions to find the quasi-stationary distribution. 
[SNl] There exists i9„ > such that 

• -0(??) < OO for < i9 < i?„, 

• ipi-d) attains its strictly negative minimum at -d* > 0, where < -d* < (and hence il''{-d*) = 
0). 

[SN2] $ is analytical in (0) for 7r/2 < < tt, where C* := iJK'^*) < 0- 

Example 13. Since for spectrally negative Levy process $(?9) is the Laplace exponent of a subordinator 
(viz. a first passage time process), the spectral representation 



(18) 



$(i9)=d_t?+/ (1 - e-''^)n_(da;), 
Jo 



applies, with J^ix A l)n_ (dx) < oo; cf. ((12). This means that if there exists a density of n_ which is 
of semixponential type, then Prop. |3] entails that $ is analytic in ^^<;* (0) and hence assumption [SN2] 
is satisfied. 

Theorem 14. If X is a spectrally negative Levy process satisfying conditions [SN1-SN2], then 

fi{dx,dy) - g_(r)2ye-''*(-+^)e"*W-K5.(x) dxdy l{,>o,,>o}, 
where (/„°° e-(*(o)+''*)^Kj* (z) dz)-i. 

Corollary 15. We have 

^lf{dy) = g_e-(*(°)+''*)-K,.(x)d.Ti{,>o} 
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and 

Observe that corresponds to an Erlang(2) distribution. The proof of these results is based on 

the following lemma, which is proven as Lemma [TOl 

Lemma 16. Under [SN1-SN2], 



asi9 I C*, where k* := y^2/ip" {■&*). 

Proof of Theorem [111 Note that all assumptions of Thm.[T]are satisfied by Prop. HI In particular, as\z\ 
tends to infinity in 5^,^* (</>) function | $(2) | is either boimded or tends to infinity. In both cases condition 
(A2) is satisfied. Moreover, for some K, 

as i9 1 C*/ ^rid we can conclude by 'Heaviside' that 
as t — > 00. Therefore, 

fl{a,l3) = lim E,[e-"Q('')"^'2W |T > t] = i_ (a)-B_ 

i— f 00 

where 



It is not hard to see that the proposed density indeed corresponds with this transform. □ 

5. Examples 

In this section we illustrate our theory by means of a number of examples. We indicate for which 
Levy processes our assumptions are fulfilled, and for a few of those processes we perform the com- 
putations. 

According to Vigon's theory of philanthropy ||T9l , a (killed) subordinator is called a philanthropist if 
its Levy measure has a decreasing density on R+ . Moreover, given any two subordinators Hi and H2 
which are philanthropists, providing that at least one of them is not killed, there exists a Levy process 
X such that Hi and H2 have the same law as the ascending and descending ladder height processes 
of X, respectively. Suppose we denote the killing rate, drift coefficient and Levy measures of Hi and 
H2 by the respective triples (6, S, 11// J and (b, 6, Il/f^). Then IT9ll shows that the Levy measure of X 
satisfies the following identity 

/•oo 

II{x,oo) — / IIho{u, oo)IIhi{x + du) + SnH-^{x) +bIlHi{x, 00), a; > 0, 
Jo 

where tthi {x) is the density corresponding to ■ By symmetry, an obvious analogue of the above 
equation holds for the negative tail n(— 00, x), with x < 0. 

Choosing then e.g. H2{t) = t and Hi{t) = T{t) with Laplace exponent $ and jump measure 11+ 
being semiexponential, then, using the above construction, we can easily give examples of spectrally 
positive Levy processes satisfying conditions [SP1-SP2]. Similarly, using the above method we can 
construct spectrally negative Levy processes satisfying [SN1-SN2]. 
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Usually these conditions can be verified in a straightforward manner, as we did in the examples 
below. 



Example 17. M/M/1 queue. In this case 
(19) 



N(t) 

X{t) = ^ a, - 
1=1 



where <Ji (where z = 1, 2, ...) are i.i.d. service times that have an exponential distribution with mean 
1/v. The arrival process is a homogeneous Poisson process N{t) with rate A; it is assumed that g := 
\/v < 1. We apply the theory of Section l4lT] 
We have 



1] + V 



= V 



7] + V 



yielding i}* — VaZ' — ly, and (* = —{\fv — \/A)^. Furthermore 

7/ + A - 1/ + ^{i-] + A - 



$(7?) = 

and hence assumptions [SP1-SP2] are satisfied. 

Let us first concentrate on We use Remark [12l Using that 

A7,2 



we obtain 



(a + d*)^'{a + rf ) - ^{a + i?*) 



This corresponds to the sum of two Exp(\/Ai') random variables. Also, 



- A 



n=0 



a + V Ai^ 



corresponding with a Geom(y^)-distributed number of Exp(v'AZ') random variables. Conclude that 
corresponds to the sum of M independent Exp(\/Ai^) random variables, where 

i.e., M has a shifted negative binomial distribution with parameters 2 and y/^. A similar form is found 
for the general light-tailed M/G/1 case. 

Let us now study /^j^^(-). It is a matter of straightforward calculus to find that 

/ . A7 - a/I \ 

i3+(/3) = (i. + /3) 



^/? + V^(V^-x/A)^ 
A partial fraction expansion argument gives that this equals 



V - \/\v + 



V - + /? 

In other words, the quasi-stationary distribution at time t (for t large) equals a mixture of an expo- 
nential and an Erlang(2) distribution. <) 
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Example 18. Linear Brownian motion. In this case X{t) = aB{t) — t, where cr > and B{t) is a standard 
Brownian motion. Remark that this process is spectrally positive and spectrally negative, so we can 
use both Thm. Eland Thm.[ll 

Let us first see what the spectrally positive results would give. It is not hard to check that 

so that, in the setting of Section |4T1 1}* = — l/cr^ and (* = —l/{2a^). It is a matter of straightforward 
computations now to obtain that 

Conclude that the quasi-stationary distributions of Q(0) and Q{t) {t large) are both Erlang(2) with 
mean whereas the stationary workload has an exponential distribution with mean 1/(20-2). 
the decomposition of Remark[l2l the first two random variables have exponential distributions with 
mean the third is equal to 0). Interestingly, the relation with the Erlang(2) distribution has also 
been observed in, e.g., 151 [131 [151 . 

The same result can be obtained by using the results from Section 14.21 Now i9* = 1 /a^ and C* = 
-l/(2cr2). It is easily checked that $(0) = 2/cr2. As expected, we obtain Jl_ (a) = A+(a) andi?_(/3) = 
B+W). 

In fact, in this case the quasi-stationarity distributions can be found in an explicit manner. For sim- 
plicity we restrict ourselves to studying just we do so by investigating the density 

:^P.(Q(0)<<z|T>t) =:ft{q). 
dq 

We rely on the standard equality 

and the fact that Q(0) (unconditioned) has an exponential distribution with mean i; here, $n( ) de- 
notes the distribution function of a standard Normal random variable, where 5'N(a;) ■— ^ — ^^lix) is 
its tail. It is known that, as a; — > oo, 

(20) vi,^(^)^fi-J_ + 3\ 1 -i.^^ 



/27r 

Let us first determine f'^(T > t) (where Q{0) has an exponential distribution with mean i), which 
can evidently be rewritten as 

Consider the first integral of the previous display. It can be evaluated as 



26-29 • dydq ^ \ / 26-2? • —^e'-^^ dqdy 

, VStT J-oo Jo v27r 

(1 - e-2(t-yVt) ) . e-^y'dy = $N(\/t) - 3'n(-\/<). 
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Likewise, the second integral can be rewritten as 



oo v27r j-oo Jo vStt 

= / (-2^-2y^/^)•^e-5y'd2/ = -2^$N(-V^) + 2^/^ 
J-oo v27r 

We arrive at 

(21) P,(r > t) = *N(\/i) + *n(-%/<) + 2i$N(-\/t) - 2Vi 

Using (|20] | it is readily verified that, for < large, 

1 -it 



/27r 



1 -it 



e 



V2^ 

Also, it holds that 

(22) ^p^(Qo<5,T>t)=2e-29vl,j,('L_^^ 2$ ^"^"'^ 



dq V / V 

so that we now have an explicit expression for ft {q), viz. the ratio of l|22t and (|2Tl l. Due to the asymp- 
totic equivalence (|20] |, Expression (|22)l behaves for t large as 

^Vf^ ^ ^ ^ — ^ 1 1 .-it 



cxp — — 7= s 2" 



^/2^ i + V 2 t y VtVi\/2^ 

We conclude that we again find that the quasi-stationary distribution of Q(0) is Erlang(2) with ex- 
pected value 2. ^ 
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